Tunneling into the edge of a compressible Quantum Hall state 
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We present a composite fermion theory of tunneling into the edge of a compressible quantum 
Hall system. The tunnel conductance is non-ohmic, due to slow relaxation of electromagnetic and 
Chern-Simons field disturbances caused by the tunneling electron. Universal results are obtained 
in the limit of a large number of channels involved in the relaxation. The tunneling exponent is 
found to be a continuous function of the filling factor v, with a a slope that is discontinuous at 
v — 1/2 in the limit of vanishing bulk resistivity p xx . When v corresponds to a principal fractional 
quantized Hall state, our results agree with the chiral Luttinger liquid theories of Wen, and Kane, 
Fisher and Polchinski. 
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Introduction: The edge of a Quantum Hall (QH) system 
plays a central role in charge transport, because the edge 
states carry Hall current Q. Also, for odd-denominator 
Landau level filling factors v that correspond to incom- 
pressible quantized Hall states, the excitations on the 
edge form a strongly interacting one-dimensional system, 
which has drawn a lot of interest [||J^]. The theoretical 
picture of the QH edge is based on chiral Luttinger liq- 
uid models, involving either one or several chiral modes 
which may travel in the same or in opposite directions. 

Another important part of the QH theory is the 
fermion-Chern-Simons approach, which can describe 
compressible QH states at even-denominator fractions 
such as v = 1/2, as well as the incompressible states j|,U. 
In this approach, the fractional QH effect is mapped onto 
the integer QH problem for new quasiparticles, compos- 
ite fermions ||, which interact with a statistical Chern- 
Simons gauge field such that each fermion carries with it 
an even number p of quanta of the Chern-Simons mag- 
netic flux. The structure of the edge can then be obtained 
from Landau levels for composite fermions in the average 
residual mag netic held @. 

The physics of the edge can be probed experimentally 
by a tunneling conductivity measurement. Soon after 
first attempts to study tunneling between edges of in- 
compressible QH states in conventional gated 2D struc- 
tures H , a new generation of 2D systems was developed 
by using the cleaved edge overgrowth technique ||. In 
these structures it is possible to study tunneling into the 
edge of a 2D electron gas from a 3D doped region. It is 
believed that the 2D electron system has a sharp edge, 
and that the the confining potential is very smooth, with 
residual roughness of an atomic scale. With the advantge 
of the high quality of the system, one can explore tunnel- 
ing into both incompressible and compressible QH states 
|jTof . It is found that the tunneling conductivity is non- 
ohmic, I ~ V a , for V > 27rfcsT/e, where the exponent a 
is a continuously decreasing function of the filling factor 



The Luttinger liquid theories predict a power law con- 
ductivity, in agreement with the experiment p0| . How- 
ever, such theories can be constructed only for specific 
incompressible densities. Furthermore, in the construc- 
tion there is no continuity in the filling factor, because 
the number of chiral modes, their propagation direction, 
and the couplings change from one incompressible state 
to the other. Thus, it is of interest to develop an alter- 
native theory capable of treating both compressible and 
incompressible states on a similar footing. 

In this paper we present a theory of tunneling into the 
QH edge based on the composite fermion picture. We 
find that under certain conditions the tunneling expo- 
nent depends only on the conductivity and interaction in 
the bulk, and is insensitive to the detailed structure of the 
edge. In this case the main effect results from the relax- 
ation of electromagnetic disturbance caused by tunneling 
electron, in which we include charge and current densi- 
ties as well as the Chern-Simons field. The characteristic 
times and spatial scales involved in the dynamics are very 
large, which makes it possible to express the leading be- 
havior in terms of measurable electromagnetic response 
functions, the longitudinal and Hall conductivities. Tun- 
neling exponents can then be found as a function of the 
filling factor. 

Let us mention that tunneling density of states in a dif- 
ferent model of a compressible edge was studied recently 
by Conti and Vignale [pOl , by Han and Thouless |0] , and 
by Han . These works focus on the case of electron 
density profile in which the density is smoothly decreas- 
ing from the bulk value iAm\k to zero at the boundary. 
It is assumed that the width of the region in which the 
density is changing is many magnetic lengths. For such 
an edge there are many hydrodynamical modes whose 
contribution to the tunneling density of states must be 
taken into account. The tunneling exponent obtained in 
this model is much larger than 1/j^uik- 
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Main results: In the calculation, we assume that the 
number of flux quanta carried by composite fermions is 
p = 2 for 1/3 < v < I and p = 4 for 1/5 < v < 1/3. Also, 
for simplicity, we assume that the composite fermions 
have "bare" conductivites pxx and pxy which are con- 
stants that may depend, for instance, on the density, but 
which arc independent of temperature. The measured 
resistivities are then p xy = p x °y +ph/e 2 and p xx — p x °x- 
The theory predicts the power-law I ~ V a , with 
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where 9jj = tan -1 p xx j p xy is the Hall angle, 6 H is the 
corresponding bare Hall angle, k is the compressibility, 
and kq — to* /2nh 2 is the bare composite- fermion com- 
pressibility, determined by the effective mass to*, which 
we treat as a constant. We have assumed, for simplic- 
ity, a short-range repulsion between the electrons (see 
Fig. 1). 
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FIG. 1. The tunneling exponent (gj) is shown as a function 
of p xy with the composite fermion flux p = 2 for 1 < p xy < 3 
and p = 4 for 3 < p xy < 5. Constant Hall angle is assumed. 
For p xx = 0.05p xy the exponent is plotted for three values 
of the model short-range interaction U = 1/k — l/«o- At 
Pxx = the exponent is universal (no {/—dependence), but at 
finite pxx it can be either bigger or smaller than the universal 
result, depending on the interaction strength kqU. 

It is interesting to compare this result with theories 
of tunneling into the edge of an incompressible quan- 
tized Hall state, where the edge is described as a single- 
or multi-channel chiral Luttinger liquid. For Jain filling 
fractions v — nj (np+ 1) with positive integer n and even 
p, there are n edge modes, all traveling in the same direc- 
tion. In this case, Wen 0] found universally a = p + 1. 
By contrast, for states such as the Jain fractions with 
negative n, where the edges have modes going in oppo- 



site directions, a may depend on the form of the interac- 
tion. Nevertheless, Kane, Fisher and Polchinski || found 
in this case that if there is sufficient scattering between 
the channels, the system scales to a universal limit, with 
a=p + l- 2/\n\. 

Our result (Q), taken at large Hall angle (Oh = 0^ = 
7r/2), agrees with the chiral Luttinger liquid theories dis- 
cussed above, when v corresponds to a Jain fraction. 
Specifically, with 



Pxy — (p + l/n)- 



Vxy ne 2 ■ 



Pxx = P W = 0, (2) 



substituted in Eq. ([!]), one gets a = 1 + \p+ l/n\ — l/|n|, 
in agreement with the universal tunneling exponents pre- 
dicted by Wen and by Kane et al. 

However, the theory presented below is not restricted 
to the incompressible states. The tunneling exponent (|l|) 
is a continuous function of the filling factor v. It is in- 
teresting, however, that a has cusp-like singularities at 
v = 1/2 and v = 1/4, in the limit p xx = 0. To under- 
stand this, consider the vicinity of v = 1/2, where the 
QH state can be described as a Fermi liquid of compos- 
ite fermions carrying two flux quanta each, and exposed 
to "residual" magnetic field 5B = 2 - v^ 1 . At v < 1/2 
the residual field direction coincides with total field, and 
all edge modes propagate in the same direction [Q. On 
the other hand, at v > 1/2, the structure of the edge is 
qualitatively different, consisting of modes going in op- 
posite directions. This effect makes v = 1/2 a singular 
density. Of course, scattering by disorder will smear the 
singularity. However, it is interesting that the change in 
the tunneling exponent ([!]) resulting from finite p xx can 
be either positive or negative, depending on the value of 
ko/ k (see Fig. 1). 

Semiclassical Green's function: We show below that 
the many-body effect on the tunneling rate arises mainly 
from the interaction of tunneling electron with slow elec- 
trodynamical modes. The low-energy QH physics can be 
described Q] by transforming electrons into composite 
fermions interacting via statistical Chern-Simons gauge 
field a^(r, t) (p = 0,1,2). At low energies, the electro- 
dynamical modes can be derived from an effective action 
S e ff[a^] for To evaluate the effect of slow fluc- 

tuations of a /J (r, t) on the electron Green's function, we 
write it as the average: 

/ D afl G rr ,(T, a^e-^ff^ 
G r * (r) = J - T (3) 



Here G rr '(ra At ) is the electron Green's function for a 
given configuration of the gauge field. For evaluating the 
tunneling current, we will only need G rr '(r) for r = r', 
which can be expressed in terms of the composite fermion 
Green's function as: 
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r(r, a fl ) = e ->° (r, a M ) 



(4) 



The exponential factor represents the action associated 
with an extra flux p&o, which must be added at t = and 
removed at t = r, to account for the difference between 
the creation operators for an electron and a composite 
fermion. Note that this factor makes the electron Green's 
function gauge invariant, while the composite fermion 
Green's function itself is not gauge invariant. 

We employ a semiclassical approximation for 
G^(t, a^). To motivate it, think of an injected elec- 
tron which rapidly binds p flux quanta and turns into a 
composite fermion. The latter moves in the gauge field 
a M and picks the phase 



■<U= / a»(r,t)jjf F \r,t)d 2 rdt, 
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i.e., at finite a xx . However, the situation is more sub- 
tle in an incompressible state with small number of edge 
modes, where there are two possibilities. First, if all edge 
modes propagate in the same direction, the orthogonality 
catastrophe effect is simply absent, because in order to 
pick up the phase Si a particle must travel across the per- 
turbation. This is the case, for instance, at Jain fractions 
v = nj (2n + 1) with n > 0, where our result is exact. On 
the other hand, if there is a small number of edge modes 
going in both directions, the orthogonality correction is 
significant and may lead to large deviations of a from 
the universal value. In this way we can understand the 
non-universal values of a obtained for Jain fractions with 
negative n in the absence of impurity scattering . 
Averaging over slow modes: Putting together Eqs. (Q)- 
(|) and Eq. § ) we get 



where (r, t) is the current describing spreading G rr (r) = G°(t) 

of free composite fermion density. Semiclassically in 
a Al (r, t) one writes 
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-S e ff[a„]+i J in a M d 2 xdt 



(7) 



G°(r) 



(6) 



where the current — r> + jf/ mJ!> , according to 



where G°(t) is the composite fermion Green's function 
in the absence of the slow gauge field. Note that fast 
fluctuations of a M are included in G (r) through renor- 
malization of Fermi-liquid parameters. Below we average 
the expression (||) over the slow gauge field (||) and find 
G rr (r) =T~ a . 

Let us discuss the physical meaning of the phase fac- 
tor approximation (^). It accounts correctly for the ac- 
tion due to charge and gauge field relaxation during tun- 
neling. One notes the similarity to the classic infrared 
catastrophe in Quantum Electrodynamics, where it has 
been shown that the phase approximation gives correct 
infrared asymptotics of electron Green's function p6[ . 
Also, similar ideas were used in the problems of tunneling 
coupled to an electromagnetic environment p7[ . 

The most important correction to (Q) is due to the 
shakeup of the composite Fermi system due to short 
range interactions with the injected electron. This "or- 
thogonality catastrophe effect" |L2| can be estimated by 
introducing scattering channels for a given field a p (r, t) . 
(For example, for the incompressible QH states the chan- 
nels are edge modes.) If the gauge field is concen- 
trated on the scale L, it excites N ~ kpL compos- 
ite fermion channels. In terms of scattering phases Si, 
(i = 1, N), the orthogonality contribution is given by 
exp(— J^i $i l 1 ^ l n (^Ao))j where to is a short time cutoff. 
To estimate Si we use Friedel sum rule Q: ^ Si = n 
(it means that we inject unit charge). Hence, Si ~ n/N 
and the orthogonality correction to the exponent a is 
~ const/iV < 1. 

The effective number of channels N involved in charge 
relaxation is certainly large in any compressible state, 
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Eq.(^), includes the gauge field flux term 

A flux \r, t) = -5, 5(r)(9(t) - 8(t r j) . (8) 



By expanding S e ff[a^] up to the second order in 
and averaging over Gaussian fluctuations in (Q), we get 
G rr (r) = G°(T)exp (S(t)) where 

S= ~ J d 3 xd 3 x'f(x)f(x') V^(x,x') . (9) 

Here V llv {x, x') = (a n (x) a v (x')) is the correlator of 
gauge field fluctuations (x = (r,t)). One may interpret 
S(t) as the semiclassical action of spreading charge [ ^8[ . 
Calculation: Hereafter we concentrate on the case of lo- 
cal conductivity tensor a, i.e., assume hi, ujI/vf < I, 
where I is the mean free path of composite fermions. 
Then the current can be found from the diffusion 

equation with a particle source at r = Jl9j 

dtjo + ViiDijVjja) = S (2 Hv)(S(t) - S(t - r)) , (10) 

Here the noninteracting composite fermion diffusion ten- 
sor D is given by the Einstein relation: D = kq^ ^ . The 
current is determined by the diffusion law j = — D Vjo- 

To get the electron Green's function at the QH state 
edge, we choose the rc-axis along the edge and assume 
that tunneling occurs at x = 0, y = (2DEG occupies 
the half-space y > 0) . Then Eq. ( [To| ) must be solved with 
the boundary condition j y (y = 0) = 0. To treat Eq. (^p|) , 
it is convenient to use Fourier transform with respect to 
x and t: jo(k,y, u>) = J jo(x, y^t) e lkx ~ luJt dxdt. One can 
write a formal solution to Eq.(|l0|): 
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J(u)S(y) , J(u) - (1 - e-™) (11) 
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Eq. ( |lip holds also in the general case if D is treated as a 
nonlocal operator. 

To find the gauge field correlator £> M „, we choose the 
gauge clq — and, following Q, write 



V 



Here 



(a'o + U~ x 
U(r - r') , 2ttU P 



(12) 



+ (13) 



is the bare interaction and Kq is current correlation 
function which can be found, e.g., using fluctuation- 
dissipation theorem: Kq = uj 2 (iuj + (VZ))V) _1 o'o . For 
simplicity, we assume short-range electron interaction 
U(r) = U5(r). 

Let us remark that one has to be careful in substituting 
( |l2| ) into (^|), since the flux current does not couple to 
the Coulomb part of the gauge field. Hence, one must 
subtract a flux term from the longitudinal current: V-j — > 

Vr* -(flux) 

After that, we decompose the product of two operators 
into partial fractions: 



'(*V)) _1 ( 
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ioj+ (V-D)V 



(14) 



iuo + (V-D)V + uj 2 & U 



(angular brackets denote integration over coordinates.) 
Now, the operators in ([u]) must be inverted, taking into 
account the boundary condition j y (y = 0) = 0. After 
some tedious algebra one finally arrives at 



S(t) = 
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1&xy k 



(15) 



where a xx and a xy are components of the measured con- 
ductivity tensor |4j] a — p^ 1 , while q 2 — k 2 + iKui/a xx , 
and <7g = k 2 + ikqoj / 'a x °J . Integrating (|l5| ) over to and k 
one gets S(t) = (a — l)ln(r/ro) with a defined by Eq. 
(|l|). Here tq is a short-time cutoff. 

We also assume that composite fermions are well de- 
fined quasiparticles and thus G® cf (t) ~ r _1 . To justify 
that, one may consider the self-energy £(e, £) of a com- 
posite fermion arising from the fluctuations of the gauge 
field. The self-energy is known to be a singular function 
of energy e but not of £ = «f(|p| — Pf) (see Q). 
Because of that, the pole of the equal point Green's func- 
tion G° r (T) remains intact. This situation is similar to 
that for electron-phonon interaction, where E does not 
depend on £ and thus does not affect tunneling current. 



We remark that the point r is taken to be close to the 
edge of the system, and the form of G° holds in this case 
even for 5B ^ 0. 

Putting everything together, we have G(t) ~ r~ a . 
The tunneling current is then given by J13|: 



I(V) - Im / G(t) 



dr - V a 



(16) 



where V is the voltage across the barrier. 

Now we check that the number of channels is big. The 
gauge field disturbance caused by the tunneling electron 
decays in the bulk as exp(— qy). Hence, after time t the 
perturbation spreads over yo(t) ~ \jK<j xx t. To estimate 
spreading along the edge, one compares two terms in 
the denominator of (15) assuming k ~ Xq , which gives 
Xo{t) — <J xy I a xx yo(t) . (Note that charge spreading is 
effectively one-dimensional if a xy 3> u xx .) To estimate 
the number of scattering channels parallel to the bound- 
ary, one takes L ~ yo{t) as the spatial scale, and gets 
A ~ yo(t)kp. Since yo(t) is large at large t (i.e., at suffi- 
ciently low bias V ~ fi/t) we indeed have N ^> 1. 
Conclusion: The above composite fermion theory of 
tunneling into a Quantum Hall edge treats on a simi- 
lar footing both compressible and incompressible states. 
We have shown that tunneling I — V characteristic is 
non-ohmic and obeys a power law for both compressible 
and incompressible Quantum Hall edges. The power law 
exponent is expressed in terms of p xx and p xy , with a 
weak dependence on the ratio k/kq, which depends on 
the form of the interaction. The exponent is a continu- 
ous function of p xy , with cusp- like sigularities at v = 1/2 
and 1/4 in the limit where p xx = 0. The physical origin 
of the singularity is that the effective magnetic field seen 
by composite fermions changes sign at v = 1/2 and 1/4. 

It is shown that the model is robust for compressible 
QH states when the number of scattering channels in- 
volved in charge relaxation is large. At the Jain fractions 
v = n/ (2n± 1) the results agree with the chiral Luttingcr 
liquid theories. 

The calculated exponent is a monotonically decreasing 
function of filling fraction, in qualitative agreement with 
the recent experiment by Chang et al jH)|,[ll]]. However, 
the theoretical value of a is larger than that observed in 
the experiments, particularly near v = 1/2. The discrep- 
ancy might be explained if the density near the edge of 
the 2D electron system is actually somewhat larger than 
that in the bulk. 
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